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Abstract 

We investigate the quantum dynamics of energy and charge transfer in a wheel-shaped artificial 
photosynthetic antenna-reaction center complex. This complex consists of six light-harvesting 
chromophores and an electron-acceptor fullerene. To describe quantum effects on a femtosecond 
time scale, we derive the set of exact non-Markovian equations for the Heisenberg operators of this 
photosynthetic complex in contact with a Gaussian heat bath. With these equations we can analyze 
the regime of strong system-bath interactions, where reorganization energies are of the order of the 
intersite exciton couplings. We show that the energy of the initially-excited antenna chromophores 
is efficiently funneled to the porphyrin-fuUerene reaction center, where a charge-separated state 
is set up in a few picoseconds, with a quantum yield of the order of 95%. In the single-exciton 
regime, with one antenna chromophorc being initially excited, we observe quantum beatings of 
energy between two resonant antenna chromophores with a decoherence time of ~ 100 fs. We also 
analyze the double-exciton regime, when two porphyrin molecules involved in the reaction center 
are initially excited. In this regime we obtain pronounced quantum oscillations of the charge on 
the fullerene molecule with a decoherence time of about 20 fs (at liquid nitrogen temperatures). 
These results show a way to directly detect quantum effects in artificial photosynthetic systems. 



* Author to whom correspondence should be addressed. Electronic mail: pulak@riken.jp 



I. INTRODUCTION 



The multistep energy-transduction process in natural photosystems begins with capturing 
sunhght photons by hght-absorbing antenna chromophores surrounding a reaction center 
[H . The antenna chromophores transfer radiation energy to the reaction center directly 
or through a series of accessory chromophores. The reaction center harnesses the excitation 
energy to create a stable charge-separated state. 

Energy transfer in natural and artificial photosynthetic structures has been an intriguing 
issue in quantum biophysics due to the conspicuous presence of long-lived quantum coher- 
ence observed with two-dimensional Fourier transform electronic spectroscopy [51 0] . These 
experimental achievements have motivated researchers to investigate the role of quantum 
coherence in very efficient energy transmission, which takes place in natural photosystems 
[5H9]. Quantum coherent effects surviving up to room temperatures have also been ob- 
served in artificial polymers |10]. Artificial photosynthetic elements, mimicking natural 
photosystems, might serve as building blocks for efficient and powerful sources of energy 
[TT| [12] . Some of these elements have been created and studied experimentally in Refs. [T51 - 
[18] . The theoretical modelling of artificial reaction centers has been recently performed in 

Refs. [niisni. 

Here we study energy transfer and charge separation in a wheel-shaped molecular complex 
(EPF complex, see Fig. 1) mimicking a natural photosynthetic system. This complex has 
been synthesized and experimentally investigated in Ref. [T7]. It has four antennas - two 
bis (phenylethynyl) anthracene (BPEA) molecules and two borondipyrromethene (BDPY) 
chromophores, as well as two zinc porphyrins (ZnPya and ZnPyb). These six light-absorbing 
chromophores are attached to a central hexaphenylbenzene core. Electrons can tunnel from 
the zinc porphyrin molecules to a fullerene F (electron acceptor). Thus, two porphyrins and 
the fullerene molecule form an artificial reaction center (ZnPya — F — ZnPyb). The BPEA 
chromophores strongly absorb around 450 nm (the blue region), while the BDPY moieties 
have good absorptions around 513 nm (green region). Porphyrins have absorption peaks 
at both red and orange wavelengths. Therefore, the BPF complex can utilize most of the 
rainbow of sunlight - from blue to red photons. It is shown in [TZj that the absorption of 
photons results in the formation of a porphyrin-fullerene charge-separated state with a life- 
time of 230 ps; in doing so, excitations from the BPEA and BDPY antenna chromophores 
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FIG. 1: (Color online) Schematic diagram of the wheel-shaped artificial antenna-reaction center 
complex reported in ref [T7]- We use the short notation, BPF Complex, to denote this photosyn- 
thetic device. The antenna-reaction center complex contains six light-harvesting pigmets: (i) two 
bis (phenylethynyl) anthracene chromophores, BPEA^ and BPEAfc, (ii) two borondipyrromethene 
chromophores, EDPY^ and BDPYf,, and (iii) two zinc tetraarylporphyrin chromophores, ZnPy^ 
and ZnPy;,. All the chromophores are attached to a rigid hexaphenyl benzene core. In addition to 
the antenna components, the photosystem contains a fullerene derivative (F) containing two pyridyl 
groups, acting as an electron acceptor. The fullerene derivative F is attached to the both ZnPy 
chromophores via the coordination of the pyridyl nitrogens with the zinc atoms. For structural 
details of the BPF Complex we refer [171 128| . 

are transferred to the porphyrins with a subsequent donation of an electron from the ex- 
cited states of the porphyrins to the fullerene moiety. This process takes a few picoseconds, 
suggesting that the excitonic coupling between chromophores is sufficiently strong. The elec- 
tronic coupling between the porphyrins and the fullerene controlling tunneling of electrons 
in the artificial reaction center also should be quite strong. It should be noted, however, 
that spectroscopic data [T514T7] show that the absorption spectrum of the BPF complex 
is approximately represented as a superposition of contributions from the individual chro- 
mophores with almost no perturbations due to the links between the chromophores. This 
means that the chromophores comprising the light-harvesting complex can be considered 
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as individual interacting units, but not as an extended single cliromophore. We can expect 
that, at these conditions, quantum coherence is able to play an important role in energy and 
charge transfer dynamics, manifesting itself in quantum beatings of chromophore popula- 
tions as well as in quantum oscillations of the charge accumulated on the fuUerene molecule. 
In principle, these oscillations could be measured by a sensitive single-electron transistor, 
thus providing a direct proof of quantum behavior in the artificial photosynthetic complex. 
Since these phenomena occurs at very short time scales (a few femtoseconds), these could 
be within the reach of femtosecond spectroscopy in the near future. The main goal of this 
study is to explore quantum features of the energy and charge transfer in a wheel-shaped 
antenna-reaction center complex at subpicosecond timescales. 

II. MODEL AND METHODS 
A. Hamiltonian 

Each chromophore has one ground and one excited state, whereas the electron acceptor 
fuUerene F has just one energy level with energy Ep. We introduce creation (annihilation) 
operators, al (a^), of an electron on the kth site. The electron population operators are 
defined as = alak- We assume that each electron state can be occupied by a single 
electron, as spin degrees of freedom are neglected. The basic Hamiltonian of the system has 
the form: 

Ho = Y^{Eknk + Ek*nk*) + EpUp + He + Vkia\,ak ajai* - ^ A^^/ala^', (1) 

k k^l era' 

where the first part incorporates the energies of the electron states (hereafter k, I = BPEA^, 
BPEA;,, BDPYa, BDPYfc, ZuPy^, ZnPyf,), and the second term is related to a fuUerene 
energy level Ep with a population operator np = Opap. The pair {k, k*) denotes a ground 
(k) and an excited (k*) state of an electron located on the site k with the corresponding 
energy Ek {Ek*). The term Hq represents the contribution of Coulomb interactions between 
electron-binding sites. This term is given in Appendix A. The fourth term of Eq. ([T]) describes 
excitonic couplings between the chromophores k and /. The matrix element Vki is a measure 
of an interchromophoric coupling strength. The last term in Eq. ([T]) describes the electron 
tunneling from excited states of the porphyrin molecules ZuPy^, ZnPyj, to the electron 
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acceptor F characterized by the tunnehng amphtudes Ao-o-', where a, a' = ZnPy*, ZnPy^, 
F. 

The interaction of the system with the environment (heat bath), represented here by a 
sum of independent oscillators with Hamiltonian 



is given by the term 



H,-p, = -T.m,u]x,,x,n,, (3) 

jk 

where Xj and Pj are the position and momentum of the jth oscillator having an effective 
mass rrij and a frequency Uj. The coefficients Xjk define the strength of the coupling between 
the electron subsystem and the environment. 

The contribution of the energy-quenching mechanisms responsible for the recombination 
processes in the system is given by the Hamiltonian 

^^quen = -J2 ill 4" «i + QiaUl* ) • (4) 

I 

For the sake of simplicity, we include the radiation damping of the excited states into the 
energy-quenching operator g/. The first term in the Hermitian Hamiltonian -ffquen is related 
to the excitation of the Z— chromophore by the quenching bath, whereas the second term 
corresponds to the reverse process, namely, to the absorption of chromophore energy by the 
bath. Both processes are necessary to provide correct conditions for the thermodynamic 
equilibrium between the system and the bath. 
The total Hamiltonian of the system is 

H = Ho + ife-ph + -f^onv + Hquen- (5) 

We omit here the Hamiltonian of the quenching (radiation) heat bath. 



B. Diagonalization of Hq 

We choose 160 basis states \M) of the complex including a vacuum state, where all 
chromophores are in the ground state and the F site is empty. We diagonalize the Hamil- 
tonian Hq ([T]) to consider the case where the excitonic coupling between chromophores, 
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described by coefficients Vim, and the porphyrin-fuUerene tunneling, which is determined 
by amphtudes A^o-', cannot be analyzed within perturbation theory. In the new basis, 
Ia*) = the Hamiltonian Hq is diagonal with the energy spectrum {E^j^}, so 

that the total Hamiltonian of the system H has the form 

H ^Y.EM{l^\ - Y.^M{^\ + ^env. (6) 

H IJ,V 

Here 

is the combined operator for both heat baths with fluctuating in time variables 

Qnu = ^mjUj'^^Xj[xjF{ti\nF\v) +Y^{xjk{n\nh\v) + (8) 

3 k 

and 

(lnv = XI (li + H.c. (9) 

I 

To distinguish the processes of energy transfer, where the number of electrons on each chro- 
mophore remains constant, from the processes of charge transfer, where the total population 
of the site changes, we introduce the following operators 

Si = ni + ni*, Ml = ni- ni*, (10) 

together with coefficients 

Thus, the environment operator Q^y can be rewritten as 

j 

with 

= J2 {^jiil^l^iW) + ^ji{l^\Mi\u)} + XjF{n\nF\i^) (13) 
I 

C. Non-Markovian equations for the system operators 

An arbitrary electron operator W can be expressed in terms of the basic operators 
Pixu = with W = Wij,^ Pfj,t,, and W^,, = The operator p^j, denotes a ma- 

trix with zero elements, with the exception of the single element at the crossing of the //—row 
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and the z/— column. The matrix elements W^u of any electron operator can be easily calcu- 
lated (see, e.g., Eqs. (SIO) and (Sll) in the Supporting Information for Ref. [SO])- For exam- 
ple, an electron localized in a two-well potential [21], with the right and left states |1) and 
|2), is described by the Pauli matrices {cxx, o'y, o"^} : = — |2)(2|, = |1)(2| + |2)(1|, 
and ay = i{\2){l\ — |1)(2|), which are expressed in terms of the basic operators |yu)(z^| with 
/i, z/ = 1,2. 

In the Heisenberg picture, the operator W evolves in time according to the equation: 
i{dW/dt) = [W,H]^ . This evolution can be described with the time-evolving operators, 
Pfj.u{t) = (1/^) which satisfy the Heisenberg equation: 



l- 



), (14) 

where = E^ — E,^, and the heat bath operator A^u is defined in Eq. ([t]). Here, we use the 
fact that the Hamiltonian H Eq. ^ is also expressed in terms of the operators p^^ taken at 
the same moment of time t. For two of these operators, p^u{t) and Pai3{t), we have simple 
multiplication rules: p^uPajs = 5i,aP^p. These rules allow to calculate commutators of basic 
operators taken at the same moment of time. We note that at the initial moment of time the 
operator, p^u{Q) = is represented by the above-mentioned zero matrix with a single 

unit at the p-z/ intersection. The matrix elements of the electron operators in Eqs. ( 9p!3 ) 
are taken over the time-independent eigenstates of the Hamiltonian Hq. The bath operators 
A^i, fluctuate in time since they depend on the environmental variables, {xj{t)}^ and on the 
variables {qiif)} of the quenching bath. 

It is known that the dissipative evolution of the two-state system can be described by the 
Heisenberg equations for the Pauli matrices {cXa,, cXy, cr^} with the spin-boson Hamiltonian 
[see Eq. (1.4) in Ref. |2T]], which includes environmental degrees of freedom. The artificial 
photosynthetic complex analyzed in the present paper has 160 states. A dissipative evolution 
of this complex is described by the Hamiltonian H in Eq. ([6]), written in terms of the 
Heisenberg operators p^y{t) = (|yu) (z/|)(t) taken at the moment of time t. Instead of the 
time-dependent Pauli matrices, the time evolution of the two-state dissipative system can be 
described by the basic operators |1)(1|, |1)(2|, |2)(1|, |2)(2|, evolving in time. In a similar 
manner, the evolution of the multi-state photosynthetic complex is described by the set 



of the time-dependent Heisenberg operators p^u{t), which obey the equation (14). As its 
spin-boson counterpart, the Hamiltonian H in Eq. ^ contains the Hamiltonian, -ffgnv, of 



the heat bath as well as the system-bath interaction terms. Here, we generalize the spin- 
boson model from the case of two states to the case of 160 states. With a knowledge of the 
operators p^jy(t), it is possible to find the time evolution of any Heisenberg operator of the 
system. Only at the initial moment of time, t = 0, the operators Pfj,u{0) form the basis of 
the Liouville space. Note that we work in the Heisenberg representation, without using the 
description based on the von Neumann equations for the density matrix. 

To obtain functions that can be measured in experiments, we have to average the operator 



p^u{t) and the equation (14) over the initial state of the electron subsystem as well as 
over the Gaussian distribution, pT = exp(— if^°^j^/T), of the equilibrium bath, (. . with 
temperature T and with a free Hamiltonian -ffbath' which is comprised of the free environment 
Hamiltonian and the free Hamiltonian of the quenching bath. The notation (...) means 
double averaging: 

(...) = ((v[/°|...|vl/0))^. (15) 

The quantum-mechanical average value of the initial basic matrix, {'^^\p^^{0)\'^'^) = 
(\&°|/i)(z/|\l/°), is determined by the product of amplitudes to find the electron subsystem at 
the initial moment of time in the eigenstates \p) and |z/) of the Hamiltonian Hq. 

A standard density matrix, p = {Pfiu}, of the electron subsystem is a deterministic 
function which allows to calculate the average value of an arbitrary operator W with the 
formula: 

{Wit)) = Tr[pit) W] =Y.W,.P^,it)- (16) 

The same average value can be written as {W{t)) = J^/nyW^^ (p^^it)), which means that 
the average matrix, {p^y{t)) = Pufj,{t), has matrix elements related to the transposed density 
matrix p{t). 

It should be emphasized that the time evolution of the heat-bath operators {xj, pj} and 
{qi}, as well as their linear combinations Qf_iu, and A^^, are determined by the total 
Hamiltonian H in Eq. ([6]). In the absence of an interaction with the dynamical system (the 
electron-binding sites), the free-phonon operators Qlfj, as well as the free operators of the 
other baths, are described by Gaussian statistics [23], as in the case of an environ- 
ment comprised of independent linear oscillators with the Hamiltonian ifenv Using the 
Gaussian property, Efremov and coauthors [23] derived non-Markovian Heisenberg-Langevin 
equations, without using perturbation theory, that assumes a weak system-bath interaction. 
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Recently, a similar non-perturbative approach has been developed by Ishizaki and Fleming 
in Ref. Due to Gaussian properties of the free bath, the total operator A^u of the 

combined dissipative environment is a linear functional of the operators p^iy, 

A,At) = + E / (^[4°iw,4p(ii)]-)^(i - ti)p-Mti). (17) 

where 9{t) is the Heaviside step function. We note that this expansion directly follows 
from the solution of the Heisenberg equations for the positions {xj} and {qi} of the bath 
oscillators. It is shown in Ref. [21] that the average value of the free operator A^^}{t) 
multiplied by an arbitrary operator B{t) is proportional to the functional derivative of the 
operator B over the variable v4^°j(t): 

{AfAmt)) =j:[dh (4°^(t)4p^(ti)) X , (18) 

fip \oAfip{ti)/ 

with 

6B{t) 



0),,, -t[B{t),pi,,{h)]^e{t-h). (19) 
Substituting Eqs. ( |l7|l8[19[ ) into Eq. (14) we derive the exact non-Markovian equation for 



the Heisenberg operators p^i, of the dynamical system (chromomorphic sites + fuUerene) 

interacting with a Gaussian heat bath, 

t 



{A'~^lit)A^^ktl)){P,ait)p,,ih)^ 

{A^^kti)Am)){p,,{t,)pUt))}. (20) 



{Pt.u)-^u;,,{p^,) = Y. j dh[{A%{tM''^lit)){pfAtl)p^^c.{t)) 



afj.u 

The time evolution of the average operator (p^;/) is determined by the second-order corre 



lation functions of the system operators as well as by the correlation functions of the free 
dissipative environment. Here we do not impose any restrictions on the spectrum of the 



environment. It should be emphasized that the exact non-Markovian equation (20) goes far 
beyond the von Neumann equation, ip = [p,H]_, for the density matrix p of the electron 
subsystem. 

D. Beyond the system-bath perturbation theory. 

We assume that the coupling of the system to the quenching heat bath determined by 
the Hamiltonian ifquen @ is weak enough to be analyzed perturbatively. However, an 



interaction of the chromophores with the protein environment cannot be treated entirely 
within perturbation theory since the reorganization energies are of the order of the intersite 
couphngs. As in the theory of modified Redfield equations EZ], the phonon operator 
Qfj,y in Eq. (12) can be represented as a sum of diagonal Q/^ = Q^J.^l and off-diagonal Q^u 
parts: 

Q^lu = Q^l5^u + (1 - 5fj,u)Q^,u- (21) 

We derive equations for diagonal and off-diagonal elements of the matrix {p^y{t)) (see Ap- 
pendix B for details about the derivation), where the interaction with the off-diagonal ele- 
ments of the environment operators Q^p are considered within perturbation theory, and the 
effects of the diagonal elements are treated exactly. 

The time dependence of the electron distribution (p^) (diagonal elements) over eigenstates 
of the Hamiltonian Hq is governed by the equation 

(Pm) + 7m (Pa.) = Il7/.a(Pa), (22) 

a 

where the relaxation matrix contains a contribution, 7^0, from the non-diagonal envi- 



ronment operators [see Eq. (B22)] as well as a contribution from the quenching processes, 



7;]^^'^ [see Eq. (B30) 



7/i« 1^la 7p 



(23) 

with the total relaxation rate 7^ = ^a'la^l■ The time evolution of the off-diagonal elements 



are given by Eq. (B31) in Appendix B. 



Equations ( 22|B31 ) allow us to determine the time evolution of an average value for an 
arbitrary operator W of the system: (yV{t)) = J2iiu{f^\^W){Pfj.u{t))- 



III. ENERGIES AND OTHER PARAMETERS 
A. Energy levels and electrochemical potentials 

The energies of the excited states of chromophores BPEA, BDPY, and ZnPy, in the 
BPF complex are estimated from an average between the longest wavelength absorption 
band and the shortest wavelength emission band of the chromophores. The average excited 
state energies of the chromophores BPEA, BDPY and ZnPy are 2610 meV, 2370 meV, and 
2030 meV, respectively, if we count from the corresponding ground energy levels [HI [T7] . 
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TABLE I: This table presents the chosen values of the excitonic couplings (V) and reorganization 
energies for energy transfer (A) of the six antenna chromophores. We choose two sets of parameters, 
one set (denoted by I) corresponds to y > A and the other set (II) to the opposite limit V < A. The 
calculated values of the time constants using both sets of parameters agree with the experimental 
values. 



Chromophores 


Set I 
Couphng (V) 


Set I 
Reorganization 
energy (A) 


Set II 
Couphng (V) 


Set II 
Reorganization 
energy (A) 


BPEAa ^ BPEAb, 
BPEA;, o BPEAa 


50 meV 


ABPEAa = 20 meV 
ABPEAfc = 20 meV 


30 meV 


Abpea„ = 40 meV 
AbpeAs = 40 meV 


BPEAa o BDPY„, 
BPEA;, ^ BDPYfo 


30 meV 


ABDPYa = 15 meV 
ARDFYf, = 15 meV 


17 meV 


Abdpy^ = 30 meV 
ABDPYi, = 30 meV 


BDPYa O ZnPya, 
BDPY;, ^ ZnPyb 


60 meV 


AznPy„ = 20 meV 
AznPyft = 20 meV 


25 meV 


AznPy, = 40 meV 
AznPy, = 40 meV 


BPEA„ o ZnPy„, 
BPEA;, o ZnPy;, 


50 meV 




40 meV 




BPEA;, o ZnPy„, 
BPEAa o ZnPyf, 


60 meV 




40 meV 





Cyclic voltammetric studies [I7j of reduction potentials with respect to the standard calomel 
electrode show that the first reduction potential of the fullerene derivative, F, is about - 0.62 
V and the first oxidation potential of ZnPy is about 0.75 V. From these data we calculate 
that the energy of the charge separated state ZnPy+— F^ is about 1370 meV. This energy 
is a sum of the energy of an electron on site F and a Coulomb interaction energy between a 
positive charge on ZnPy and a negative charge on F. The Coulomb energy can be calculated 
with the formula u = e^/47reoer, where eo is the vacuum dielectric constant. The dielectric 
constant e of 1,2 difiurobenzene (a solvent used in all experimental measurements of Ref. |17j ) 
is about 13.8. If the distance r between porphyrin ZnPy and fullerene F is about 1 nm, the 
Coulomb interaction energy is about 105 meV. Thus, the estimated energy of the electron 
on F can be of the order of 1475 meV. 
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B. Reorganization energies and coupling strengths 

The reorganization energies for exciton and electron transfer processes and electronic 
coupling strengths between the chromophores depend on the mutual distances and orienta- 
tions of the components, strengths of chemical bonds, solvent polarity and other structural 
details of the system. Precise values of these parameters are not available. However, time 
constants for energy transfer between different chromophores in the EPF complex, as well as 
rates for transitions of electrons between the fuUerene F and porphyrin chromophores ZnPy, 
have been reported in Ref. jT^. We fit the experimental values of these time constants with 
the rates following from our equations with the goal of extracting reasonable values for the 
reorganization energies and the electronic and excitonic couplings. In principle, many com- 
binations of reorganization energies and coupling constants could be possible. For the sake 
of simplicity, we consider two sets of parameters, for two limiting situations. One parameter 
set (denoted by I in Table I) corresponds to a larger excitonic couplings, V, compared to 
the reorganization energies, A, whereas another set of parameters (denoted by II in Table I) 
considers the opposite case: where the reorganization energies are larger than the excitonic 
couplings. These two sets of parameters are presented in Table I. In addition to the parame- 
ters listed in Table I, we consider the following values for the charge-transfer reorganization 
energies (set I): Xp = 200 meV, Xim = 100 meV, and Xp = 230 meV, Xim = 120 meV (set 
II), where / = ZnPy^, ZnPyj,. The values of the reorganization energies for energy-transfer 
processes are much smaller than those for charge transfer. 

References [171 |28] reported a very fast electron transfer (with a time constant r ~ 3 
ps) between excited states of zincporphyrins (ZnPya,ZnPyb) and the fuUerene derivative F. 
This fact indicates a good porphyrin-fuUerene electronic coupling, which is due to the short 
covalent linkage and close spatial arrangement of the components [28j . Hereafter, we assume 
that the ZnPy-F tunneling amplitudes A are about 100 meV (parameter set I) and 80 meV 
(parameter set II). These parameters provide a quite fast electron transfer, despite of a 
significant energy gap between the ZnPy excited states and the fuUerene energy level. 

To describe recombination processes, we introduce a coupling of the l-th chromophore 
to a quenching heat-bath characterized for simplicity by the Ohmic spectral density: 
Xii^) =0-1^ with a dimensionless constant a/. We assume that the shifts of the en- 
ergy levels caused by the quenching bath are included into the renormalized parameters 
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TABLE II: This table presents a comparison between the calculated values of the time constants 
(using the parameters sets I and II) to the experimental values reported in Ref. [17]. 



Process 


r (Set I) 


r (Set II) 


r (Experimental) 


BPEAa BPEAb, 
BPEAft ^ BPEAa 


~ 0.4 ps 


~ 0.2 ps 


0.4 ps 


BPEA„ BDPY„, 
BPEAfe ^ BDPYb 


~ 5 ps 


~ 5.4 ps 


5-13 ps 


BDPY„ ZnPv«, 
BDPYfc ZnPyb 


~ 5 ps 


~ 3.9 ps 


2-15 ps 


BPEA„ ZnPvo, 
BPEAfo ^ ZnPyb 


~ 12 ps 


~ 12 ps 


7 ps 


BPEAb ^ ZnPy,, 
BPEA„ ^ ZnPyb 


~ 10 ps 


~ 12 ps 


6 ps 


ZnPyfe ^ F, 
ZnPy, ^ F 


~ 3 ps 


~ 3 ps 


3 ps 



of the electron subsystem. The experimental values [T71 [28] of the lifetimes rf for excited 
states of chromophores BPEA, BDPY and ZnPy: r^p^A = 2.82 ns, T^ppY = 0.26 ns, 
and t|jjpy = 0.45 ns, can be achieved with the following set of coupling constants: 
"BPEA ~ 10"^, aBDPY ~ 10~^ and aznPy ~ 7 x lO"'^. 



IV. RESULTS AND DISCUSSIONS 



and two sets of parameters discussed in Sec. Ill, here we study 
electron and energy transfer kinetics in the BPF complex with special emphasis on the 
femtosecond time range, where the effects of quantum coherence can play an important role. 
We consider both single- and double-exciton regimes. 



A. Evolution of a single exciton in the BPF complex 

In Fig. 2 we show the time evolution of the excited states populations provided that only 
the BPEAq chromophore is excited at t = (single-exciton regime). We use here the param- 
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Using Eqs. (B3Tp2 
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FIG. 2: (Color online) Site populations as a function of time for the parameter set I. The inset 
plots depict the features of site populations for short times, at two different temperatures: T = 
300 K and 77 K. The site populations of the BPEA moieties oscillate with a considerably large 
amplitude, while the oscillations of the other site populations are hardly observable. 



eter set I, where excitonic couplings are larger than reorganization energies (see Sec. III). 
The process starts with quantum beatings between the resonant BPEAq and BPEAf, chro- 
mophores, with a decoherence time of the order of 100 fs (at T = 300 K). In a few picoseconds, 
the excitation energy is subsequently transferred to the adjacent BDPY moieties and to the 
ZnPy chromophores. Later on, an electron moves from the excited energy level of the por- 
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FIG. 3: (Color online) This figure presents site populations as a function of time for the parameter 
set II. The inset plots show the site populations for short times, at two different temperatures: T 
= 300 K and 77 K. The amplitudes of the site-population oscillations are much smaller and die out 
earlier, compared to Fig. 2. This figure indicates that even for A > y, the energy transfer between 
BPEA chromophores is dominated by wave-like coherent motion. 

phyrins to the fullerene moiety; thus, producing a charge-separated state, ZnPy+— F~, with 
a quantum yield 95%, which is in agreement with experimental results [I6]. It is evident 
from Fig. 2 that excited state populations of the BDPY chromophores oscillate with much 
lower amplitudes and die out within a very short time, t < 10 fs, at both temperatures: T 
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FIG. 4: (Color online) Site populations as a function of time for the parameter set I, when the 
ZnPya chromophore is in the excited state and all the other chromophores are in the ground state 
at t = 0. The inset plots depict the site populations at short times for two temperatures: T = 300 
and 77. Lowering the temperature enhances the oscillations of the charge density on the fullerene 
moiety. Despite the huge energy difference between ZnPy*— F and ZnPy+— F~, the charge of the 
fullerene site exhibits oscillatory behavior for short times, specially at lower temperatures. 

= 300 K and 77 K. The populations of the other sites of the BPF complex do not exhibit 
any oscillatory behavior. This can be ascribed to incoherent hopping becoming dominant 
because of significant energy mismatch between these chromophores. 

Figure 3 shows the time-dependence of the excited state populations of chromophores for 
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the parameter set II, where the reorganization energies are larger than the excitonic couplings 
between chromophores. At t = the BPEAq chromophore is excited (single-exciton regime). 
Then, after a few picoseconds, the charge-separated state is formed with a quantum yield 
of the order of 97%. However, owing to a stronger system-environment coupling, quantum 
beats between the BPEAq and BPEAf, chromophores have a lower amplitude and shorter 
decoherence time (~50 fs) than in the previous case when we use the parameter set I. We 
note that no quantum oscillations of the fuUerene population (site F) are visible in Figs. 2 
and 3. 

No significant oscillations of the site populations were observed (not shown here) when 
the BDPY chromophores were initially (at t = 0) excited. In this case, due to the consid- 
erable energy gaps between the BDPY and the adjacent BPEA and ZnPy chromophores, 
incoherent hopping dominates over the coherent transfer of excitons. Furthermore, the struc- 
ture of the BPF complex [151 12H] does not allow direct energy transfer between two BDPY 
chromophores. 

Figures 4 and 5 demonstrate charge- and energy-transfer dynamics for two parameter sets, 
I and II, for the case when one of the porphyrin chromophores (ZnPy^) is excited. Here we 
do not show the time evolution of the BPEA and BDPY chromophores since these moieties 
have higher excitation energies than the ZnPy chromophore and they are not excited in the 
process. As evident from Figs. 4 and 5, the excited porphyrin molecule rapidly transfers an 
electron to fullerene, thus, producing a charge-separated state ZnPy+— F^ with a quantum 
yield of about 98%. The most important feature here is that the population and charge of 
the fullerene molecule oscillates in time due to a quantum superposition of the porphyrin 
excited state and the state of an electron on the fullerene. The amplitude of these quantum 
beats is very small and the decoherence time is quite short (~10 fs at T = 77 K). This fact 
can be explained by the significant energy mismatch between the ZnPy*— F and ZnPy"^— F~ 
states as well as by the strong influence of the environment on the electron dynamics. 

B. Evolution of double excitons in the BPF complex 

In the previous subsection, we consider a single exciton case with just one chromophore 
initially being in the upper energy state. Here we analyze a situation where two porphyrin 
molecules (ZnPy^ and ZnPyf,) are excited at t = 0. Figures 6a and 6b show the coherent 
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FIG. 5: (Color online) Time evolution of the site populations for the parameter set II, starting 
with an exciton on the chromophore ZnPya at t = 0. The inset plots depict the features of the site 
populations for a shorter time regime and at two temperatures: T = 300 K and 77 K. Lowering 
the temperature enhances oscillations of the charge density on the fullerene derivative. These 
results indicate that the population of the site F oscillates for short times, even for A > V. These 
oscillations are more pronounced at lower temperatures. 

dynamics of the fullerene population (and the fullerene charge) for the parameter sets I 
(Fig. 6a) and II (Fig. 6b) at two different temperatures, T = 77 K and T = 300 K. We 
also compare the double-exciton case with the previously analyzed single-exciton case. It 
is apparent from Fig. 6, that the double excitation significantly enhances the amplitude 
of quantum oscillations of the fullerene charge for both sets of parameters. As one might 
expect, the frequency of the quantum beatings and the decoherence time are not affected 
by the number of excitons. 
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FIG. 6: (Color online) Time evolution of the populations on the site F, for both sets of parameters, 
I and II, comparing the double-exciton case (the two ZnPy chromophores are excited) with the 
single-exciton case, (a) Time evolution of the populations on the site F for the parameter set I. 
(b) Time evolution of the populations on the site F for the parameter set II. Note that the double- 
excitation significantly enhances the amplitude of the charge oscillations at the fullerene site for 
both sets of parameters, either at low or high temperatures. 

C. Amplification of charge oscillations 



In the previous discussion we observed that lowering the temperature and the simultane- 
ous excitation of both porphyrins significantly enhances quantum oscillations of the fullerene 
charge. In this subsection we show that these oscillations can also be controlled by tuning 
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Time (fs) 

FIG. 7: (Color online) Time evolution of the population on the site F for the parameters set II 
when both ZnPy chromophores are excited at t = 0. (a) Effects of the coupling A on the time 
evolution of the populations on the site F. (b) Effects of the energy gap between an excited state 
of a ZnPy chromophore and the charge-separated state, E^h, on the time evolution of populations 
on the site F. (c) Effects of the reorganization energy A on the time evolution of populations on the 
site F. As can be seen from these plots, the contribution of wave- like coherent motion to electron- 
transfer dynamics is significantly enhanced when strengthening the coupling between fullerene and 
porphyrin, lowering the energy gap between the fullerene and porphyrin sites, and decreasing the 
reorganization energy. 
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the following parameters: 

1. Electron tunneling amplitude A . 

The electronic coupling between the fullerene electron acceptor and zinc porphyrins has 
a strong effect on the quantum oscillations of the fullerene charge. To explore this effect, in 
Fig. 7a we plot the electron population of the fullerene as a function of time, for different 
values of the coupling A. Figure 7a clearly shows that, with increasing A, the amplitude of 
the charge oscillations is significantly enhanced. This coupling can be increased by attaching 
the fullerene to porphyrins with better ligands which form much stronger covalent bonds. 

2. Energy of the charge-separated state Ech ■ 

The energy Ech ~ 1370 meV, of the charge separated state, ZnPy+— F^ is much lower 
than the energy of the zinc porphyrin excited state, EznPy* ~ 2030 meV. It is evident from 
Fig. 7b that increasing the energy E^h, which leads to a decrease of the porpyrin-fuUerene 
energy mismatch, results in a pronounced amplification of the quantum oscillations of the 
fullerene charge. The energy of the fullerene can be changed by placing nearby a charge 
residue, electrostatically coupled to the fullerene. 

3. Reorganization energy Xp . 

In Fig. 7c we present the time evolution of the fullerene population for different values of 
charge transfer reorganization energy Xp. This parameter can be decreased by replacing the 
polar solvent with another one which has a much lower polarity. As can be seen from Fig. 7c, 
the quantum oscillations of the fullerene charge survive much longer times for smaller values 
of the reorganization energy, which correspond to weaker system-environment couplings. A 
similar effect is expected when the porphyrin reorganization energy is changed. 

V. CONCLUSIONS 

We theoretically studied the energy and electron-transfer dynamics in a wheel-shaped 
artificial antenna-reaction center complex. This complex [T7], mimicking a natural photo- 
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system, contains six chromophores (BPEAa, BPEAf,, BDPYa, BDPYfc, ZnPy^, ZnPyt) and 
an electron acceptor (fullerene, F). Using methods of dissipative quantum mechanics we 
derive and solve a set of equations for both the diagonal and off-diagonal elements of the 
density matrix, which describe quantum coherent effects in energy and charge transfer. We 
consider two sets of parameters, one corresponding to the case where the energy-transfer 
reorganization energy A is less than the resonant coupling V between the chromophores, 
A < l^, and another regime where A > V. For these two sets of parameters we examine the 
electron and exciton dynamics, with special emphasis on the short-time regime (~ femtosec- 
onds). We demonstrate that, in agreement with experiments performed in Ref. [T7], the 
excitation energy of the BPEA antenna chromophores is efficiently funneled to porphyrins 
(ZnPy). The excited ZnPy molecules rapidly donate an electron to the fullerene electron 
acceptor, thus creating a charge-separated state, ZnPy+— F^, with a quantum yield of the 
order of 95%. There is no observable difference in energy transduction efficiency for these 
two sets of parameters. In the limit of strong interchromophoric coupling, coherent dynam- 
ics dominates over incoherent-hopping motion. In the single-exciton regime, when one of 
the BPEA chromophores is initially excited, quantum beatings between two resonant BPEA 
chromophores occur with decoherence times of the order of 100 fs. However, here the elec- 
tron transfer process is dominated by incoherent hopping. For the case where one porphyrin 
molecule is excited at the beginning, we obtain small quantum oscillations of the fullerene 
charge characterized by a short decay time scale (~ 10 fs). More pronounced quantum os- 
cillations of the fullerene charge (with an amplitude ~ 0.1 electron charge and decoherence 
time of about 20 fs at T = 77 K) are predicted for the double-exciton regime, when both 
porphyrin molecules are initially excited. We also show that the contribution of wave-like 
coherent motion to electron-transfer dynamics could be enhanced by lowering the temper- 
ature, strengthening the fuUerene-porphyrin bonds, shrinking the energy gap between the 
zinc porphyrin and fullerene moieties (e.g., by attaching a charged residue to the fullerene), 
as well as by decreasing the reorganization energy (by tuning the solvent polarity). 
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Appendix A: Coulomb interaction energies 

The Coulomb interactions between the electron states are, 

He = -Uf [(1 - nznPyJnp + (1 - nZnPyb)^F] + MPy(l " nZnPyJ(l - ^ZnPyb) 

+ MZnPya'^ZnPya'^ZnPyJ + MZnPyt'^ZnPyb'^ZnPy*, (Al) 

where, 

^ZnPya = ""-ZnPya + ^ZnPyJ and ?^ZnPyb = ""-ZnPyb + ^ZnPyf • 



The first term of (Al) represents the electrostatic attraction (so the minus sign) between 
the positively charged ZnPy chromophores and the negatively-charged fuUerene. The sec- 
ond term is due to the Coulomb repulsion (so the plus sign) between two ZnPy chro- 
mophores. The last two terms are the repulsive interaction energies when both the excited 
and ground states of the ZnPy chromophores are occupied by electrons. The coefficients 
mf, Mpy, MznPya; cLud MznPya represent the magnitude of the electrostatic interactions and 
these are calculated using the Coulomb formula. We have assumed that the empty ZnPy 
chromophores (nznPy + '^znPy* = 0) have positive charges and the acceptor state F becomes 
negatively-charged when it is occupied by an electron. 

Appendix B: Derivation of equations for the matrix (p^u) 

Our derivation of the equations for the matrix (p^u) is based on the exact solution for 
the operator p^^ = (|/i)(z/|)(t) of the system influenced only by diagonal fluctuations of the 
bath. In this case the "system -|- bath" Hamiltonian has the form 



where = A^^ [see Eq. (13)]. The time evolution of the exciton operators is governed 



by the Heisenberg equation 

^P^lu = - i^^iuPpiu + E ^j^'ji^j - ^'j)xjP^,u■ (B2) 
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It is possible to verify that the solution of Eq. (B2) is given by the equation 



Piiu{t) = exp[iQ^^{t - to)] X exp 



exp 



P^lu(to), 



where 



E 



u\2 



(B3) 



(B4) 



and pj is the Heisenberg operator of the dissipative environment. The evolution begins at 
time t = to. The diagonal operators = p^^ are constant, Pfj,{t) = p^{to), in the presence 
of a strong interaction with the diagonal operators of the protein environment. 

For uncorrelated diagonal and off-diagonal environment operators, when 
{Qa\^)Q^^K^')) ~ 0) contribution of the environment to the non-Markovian equation 



(20) consists of two parts: 



(B5) 



The diagonal elements, Q^, of the environment contribute to the first part, 

•J 

dU{Qf{U){Qf -Qf){t)){p,{U)p,^{t)), (B6) 



whereas the second part is due to a contribution of the non-diagonal (abbreviated as n-diag 
in the super-index) operators, Q^y, 



dti{&)Qf^ih)){p,^it)pp,ih)) + 



dti{Qfhti)Ql-I{t)){p;i,{ti)p,a{t)) + 

dh{Q^^l{t)Q^^kti)){Pa.{t)pMti)) - 
dh{Q^^kti)Q'^'lit)){pMti)po.At)). (B7) 

We note that the time evolution of the diagonal elements of the system operator, p^ = p^^, 
is determined by the non-diagonal operators Q^i^ as well as by quenching terms. Strong 
diagonal fluctuations of the environment have no effect on the evolution of the diagonal 
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elements of the matrix. Thus, in Eq. (B6) we assume that pp(ti) = Pp(^), so that Eq. (B6) 
can be rewritten as 



(B8) 



where the time-dependent rate, r^'^s(t), and the frequency shift, 5^2^'^^, can be found from 
the following expression 

J 

The rate '^^^^^if) determines the fast decay of quantum coherence in our system. For an 
environment composed of independent oscillators we obtain 

- E ^(A.^ - A coth (^) cosa;,(t - t,) - 



(A^)2-(An' sina;,(t-t 



(BIO) 



The fluctuations of the diagonal operators of the environment can be described by the set 
of spectral functions, 

,3 



(Bll) 



together with the corresponding reorganization energies. 



A 



(B12) 



We also introduce a spectral function, J^^^co), which characterizes the non-diagonal (/i 7^ u) 
environment fluctuations, 

^..M = E ^|ArP5(a; - c,), (B13) 



where = h.^^ (13) taken at /i 7^ z/. With Eq. (BIO) we calculate the contributions of the 
diagonal environment fluctuations into the decoherence rate and the frequency shift of the 
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off-diagonal elements of the system matrix {pfj^,) in (B8) 



J^,y{u}) coth ( ) smut, 
u \2T J 



tu 



(B14) 



The contribution of the non-diagonal fluctuations of the environment to the evolution of 



the electron operators {Pfj.^) is defined by Eq. (B7). To calculate the products of exciton 



variables taken at different moments of time, for example, p^aif)pjiuiti)^ we use Eq. (B3), 
which describes the evolution of exciton operators in the presence of strong coupling to the 
diagonal operators, Q^, of the environment. We assume that the interaction with the non- 
diagonal environment operators, Q^v, is weak. With Eq. (B3) we express the operators at 
time ti in terms of operators taken at time t: 



Pfiuih) = exp [-iQfipT] exp [iupp^r)] exp [-w^p(t, h)] p,_,u{t), 
Pppih) = Pi,u{t) exp [-iVtjii^T] exp [-m^p(r)] exp ti)] 

where t = t — ti, and 



(B15) 



duj 



w 



■J^^(w) sin wr. 



(t,ti) = E(A'-Ai)[P.W-P.(^i)]- 



(B16) 



Here we assume that pj{t),pj{ti) are free-evolving momentum operators of the environment, 
which are described by Gaussian statistics with a correlation function 
1 



hrrijUj , / hojj \ , , 

coth — ^ cos Ui(t — ti). 

2 \2T 



(B17) 



The operator function Vfj,u(t,ti) does not commute with the exciton matrix p^uit), and, 
therefore, we need two expressions for the operator p^uiti), which are distinguished by the 
order of the operators p^uif) and exp [— it>^,^(t, ti)] . For the average value of the operator 
exp [—iv^i,{t,ti)\ we obtain 

(exp = exp ^J^^(u;) coth [1 - cosa;(t - . (B18) 



Substituting Eqs. (B15) to Eq. (B7) and using the secular approximation we obtain a 



contribution of the non-diagonal environment operators, Q^j,ui to the evolution of diagonal 
exciton operators (p^), 

(-^[P.,^e%'r]-) = -E7..(t)(p.) +E7m«W(p«), (B19) 
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characterized by the following relaxation matrix, 



) + 



-tl) iUafj.{t-tl) I -iVaf_i{t,t-L)\ 



(B20) 



where 



= (1/2) / Ja,{u^) X 



coth ( ^] - 1 



Juj(t-ti) _|_ 



coth(|,)+l 



(B21) 



When the environment is at high temperatures (2T ^ cj) and at low frequencies of the 
diagonal fluctuations {ojt <^ 1) we have: 

and 

{exp[-iv^u(t,ti)]) ~ exp[-A^,,r(t - tif]. 
With these assumptions the relaxation matrix has a simple form 



7r 



Jo 



duj Jafii^^) n{uj) X 



exp 



exp I - ) exp 



4A,^r 



(B22) 



where 71(0;) = [exp((X'/T) — 1] ^ is the Rose distribution function at the temperature T. 



The moment of time t in the expression (B20) for the relaxation matrix is usually higher 



than the effective retardation time, Tc ~ (Xa^T) of the integrand in Eq. (B20): t ^ Tc. 
Therefore, we assume that t ^ 00, so that ^^a(t) — 7/^q(oo) = j^a- 



It follows from Eq. (B7) that a contribution of the non-diagonal environment operators 



Q^i, to the evolution of the off-diagonal elements p^i, is given by the formula 



(B23) 



where 
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). (B24) 



A small frequency shift, Sfln^, can be hereafter ignored. The dephasing rate, f ^j,, has two 



parts, = + Tj,, where 



TT 



Xfj,aT Jo 



exp 



CO 



+ exp ( - ) exp 



(B25) 



We note that = (1/2) Ea 7^^, and fi^^ = + from Eq. (|B4[),(|B12[). 

Assuming that the environment fluctuations acting on each electron-binding site are 



independent and using Eq. (13) for the coefficients A^'', we obtain 



J^,{uj) = [jis{uj)MSi\u)\^ + Jim{uj)MMi\u)\^] + JF{uj)MnF\u)\ 



(B26) 



where 



Jf(uj) 



(B27) 



The results obtained above are valid for an arbitrary frequency dependence of the spectral 
densities Jis{uj), Jim{uj), Jf{oj). Hereafter we assume that these functions are described by 
the Lorentz-Drude formula characterized by a common inverse correlation time, 7c = t~^, 
and by a corresponding reorganization energy \is, Xim, or \f, e.g. 

^ Ms W7c 



(B28) 



Quenching processes also contribute to the decay of the off-diagonal elements, {p^j^u), with 



the following decoherence rates: T^'^^ 



■pquen 



la 



COth ( ^ 

2T 



+ 1 



(B29) 



Here we consider an Ohmic quenching heat-bath with the spectral density x'/(w) = aiu, 
which is determined by a set of site-dependent dimensionless coupling constants a; ^ 1. 
The contribution of quenching to the relaxation of the diagonal elements of the electron 
matrix, (p^), is determined by the standard Redfield term 



7^^^ = E(K/^l4aH^)r + mala,,\p)\')x':{u 



flU J 



cothl^)-l 



(B30) 
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As a result, we find that the time evolution of the off-diagonal elements of the electron 
matrix is determined by the expression 



{pi^u){t) = exp{iujf,^t- X^^Tt^) X exp(-r^^t) p^y{Q), 



(B31) 



with the decoherence rates T uy = r„ + Ti,, where the coefficient r„ contains contributions of 



the off-diagonal fluctuations of the environment (B25) as well as quenching processes r^"'^'^ 
(B29): = + r^"^'^. The evolution starts at the moment t = with the initial matrix 



Pij.iy{0)- An effect of diagonal environment fluctuations is determined by the rate JX^^T 



where A^j, is the reorganization energy defined by Eq. (B12) and T is the temperature of 
the environment. 
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